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INTRODUCTION 

Reliable estimates of stability and control derivatives are essential 
for flight simulations and handling quality evaluations of aircraft. In 
response to the growing need for reliable derivative estimates, the NASA 
Dryden Flight Research Center developed a technique for obtaining the sta- 
bility and control derivatives of aircraft from flight data (ref. 1) and 
developed a set of FORTRAN computer programs to implement the technique 
(ref. 2). This method of derivative extraction is based on a modified 
maximum likelihood estimator that uses the Newton-Bal akrishna'*. algorithm 
to perform the required minimization. 

These computer programs were used to determine the stability and control 
derivatives of a modified Beech 99 airplane. The aircraft, flown as a co- 
operative effort by NASA, Beech Aircraft Corporation, and the University 
of Kansas Flight Research Laboratory, was utilized to study the effects 
of separate surface stability augmentation (ref. 3). Data were obtained 
with the aircraft in a clean configuration and with one-third flap deflec- 
tion. This report presents the Beech 99 derivative estimates obtained 
with the modified maximum likelihood technique and compares these esti- 
mates with predicted values. 


SYMBOLS 


A, B, C, 0, R 



System matrices 

normal, longitudinal, and lateral 
accelerations, g 

reference span, m 

coefficients of lift and drag 

coefficients of roll, pitch, and 
yaw moment 

coefficients of normal and lateral 
force 

center of gravity 

mean aerodynamic chord, m 

measurement noise spectral 
density matrix 

acceleration due to gravity, m/sec 2 
moments of inertia, kg - m 2 
cost functional 
mass, kg 

roll rate, deg/sec 
pitch rate, deg/sec 
dynamic pressure, N/m 2 
yaw rate, deg/sec 
wing area, m 2 
time, sec 

thrust coefficient 
total time, sec 
control vector 
velocity, m/sec 



Superscripts: 


weight, kN 


X, V, 2 


longitudinal, lateral, and normal 


XAN, XAY 


distances of an and ay accelerometers 
forward of center of gravity, m 


state vector 


computed observation vector 


distance of a^ accelerometer 
below center of gravity, m 


measured observation vector 
angle of attack, deg 
angle of sideslip, deg 


6 a’ V 6 r 


aileron, elevator, and rudder 
deflections, deg 


Subscripts: 


measurement noise vector 


pitch attitude, deg 
vector of unknowns 
roll attitude, deg 


p, q, r, a, 0, 6 a , 6 e , 


derivative with respect to indicated 
quanti ty 


matrix transpose 
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DESCRIPTION OF THE AIRPLANE AND INSTRUMENTATION 


The modified Beech 99 airplane used in this analysis is a 14-feat, twin 
turboprop commercial airliner with low wings and retractable landing gear 
(figs. 1 and 2). Tables 1 and 2 list important geometric and mass charac- 
teristics of the Beech 99 airplane. The test airplane was modified so that 
there were two independently operable control surfaces where there is normal- 
ly only one; however, only one set of rudder, aileron and elevator surfaces 
was used during the test flight. The extra surfaces are shaded in figure 1. 
These extra surfaces remained in fixed positions during the flight. 

The instrumentation of the airplane consisted of a standard package 
used for the measurement of stability and control parameters, including 
three-axis angular rate gyros, attitude gyros, and linear accelerometers, 
along with boom-mounted angle of attack and angle of sideslip es. The 
data were filtered with 40-hertz passive analog filters, then stapled with 
a 9- bit pulse code modulation (PCM) system and telemetered to a ground 
station for real-time monitoring. 

The analysis used in the derivative extraction accounts for the effect 
of instrument location on the measurement of linear accelerations and flow 
angles. The instrument locations used in the analysis of the flight data 
are presented in table 3. Table 4 lists the resolutions of the instrumen- 
tation system used in the analysis. 


TEST PROCEDURES AND FLIGHT CONDITIONS 


The Beech 99 airplane was flown in the cruise configuration for half of 
the maneuvers analyzed and at one-third flap setting for the remainder. All 
maneuvers were flown with the center of gravity at approximately 26 percent 
of the mean aerodynamic chord. Most of the maneuvers performed were simple 
aileron, rudder, or elevator pulses. 

All data were obtained during a 2-hour flight in smooth air. Fifty-six 
maneuvers were performed for derivative estimation over an angle of attack 
r w e from 1.5 to 4.5 degrees, a velocity range from 65 to 100 meters per 
second and an altitude range from 1800 to 3200 meters. Table 5 lists the 
flight condition corresponding to each maneuver. 


METHOD OF ANALYSIS 


A maximum likelihood method of analysis was used to determine a complete 
set of linear body axis stability and control derivatives from the 56 maneu- 
vers performed in flight. The digital computer program used is called the 
modified maximum likelihood estimator, version three (MMLE-3) which is an out- 
growth of the program described in detail and listed in reference 2. The 
program is discussed briefly in appendix A.. Further information is available 
upon request from the Dryden Flight Research Center. The analysis technique 
is an iterative technique that minimizes the difference between the measured 
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aircraft response and the computed aircraft response by adjusting th ■ stabil- 
ity and control derivative values used in calculating the computed r<* rinse. 
This method can be modified to include a priori information from prev'.rs 
calculations, flight tests, or wind tunnel tests; however, no a priori infor- 
mation was used in this Beech 99 analysis. The maximum likelihood techni- 
que is described fully in reference 1. 

In addition to giving derivative estimates, this method provides uncer- 
tainty levels for each derivative. The uncertainty levels are proportional 
to the approximation of the Cramer-Rao bounds described in reference 1, and 
are analogous to the standard deviations of the estimated derivatives. The 
larger the uncertainty level, the more uncertain the validity of the estimated 
value. The uncertainty levels obtained for derivatives from different maneu- 
vers at the same flight conditions can be compared to determine the most valid 
estimate. The uncertainty levels provide additional information about the 
validity of the derivative estimation. Further information on the interpre- 
tation of uncertainty levels is included in reference 4. 

The digital computer program used in the data analysis is capable of pro- 
ducing one set of derivative estimates based on multiple sets of data. The 
lateral-directional derivative estimates in this report result from the simul- 
taneous analysis of both rudder and aileron maneuvers. Analysis of rudder 
and aileron maneuvers simultaneously usually results in improved derivative 
estimates, as shown in reference 4. 

The cost functional minimized by the program is given in appendix A. 

The matrix G in this cost functional acts as a signal weighting matrix. For 
this analysis , G was chosen to be diagonal with values given in table 6. 


RESULTS AND DISCUSSION 


For all 56 maneuvers flown with the Beech 99 airplane, the measured air- 
craft response compared satisfactorily with the computed response based on 
the maximum likelihood estimation. A typical longitudinal maneuver is shown 
in figure 3, and a typical lateral-directional double maneuver is shown in 
figure 4. The measured (solid-line) and computed (dashed line) response of 
the aircraft are in excellent agreement in these figures. Some maneuvers pro- 
duced better agreement than others; however, on the average, the agreement 
was quite good. 

An atypical maneuver combining aileron and rudder inputs is shown in 
figure 5. This maneuver was performed to determine any nonlinearities in 
aileron effectiveness. The excellent fit in figure 5 based on a linear model 
indicates that there was no significant nonlinearity in aileron effectiveness. 
The MMLE stability and control derivative estimates based on the analysis of 
this maneuver were in excellent agreement with estimates from more conven- 
tional maneuvers. The analysis of this maneuver required that the small angle 
approximation (in bank angle) be removed from the mathematical model (see 
appendix B). All other maneuvers were analyzed using small angle approxima- 
tions. 

The estimates for 56 maneuvers (19 primarily with elevator inputs, 21 
primarily with rudder inputs, 15 primarily with aileron inputs, and one with 
simultaneous aileron and rudder inputs) are presented in figures 6 and 7. 
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The longitudinal stability and control derivative estimates are pres- 
ented in figure 6, while the lateral -directional estimates are displayed in 
figure 7. Each symbol indicates the derivative estimate for one maneuver, 
and the vertical bar associated with the symbol represents the uncertainty 
level for that derivative estimate. The square symbols are used to indicate 
one-third flap down maneuvers. Manufacturer predicted derivatives based on 
a combination of analytical predictions, wind tunnel tests, and flight test 
results are shown along with the MMLE-3 derivative estimates for comparison. 

The manufacturer’s predictions of derivatives that are functions of 
thrust coefficient (C N , C N , C m , and C n ) are indicated by dashed lines 

« 6 e 6 e P 

for thrust coefficients of 0 and 0.1. This range of thrust coefficients in- 
cludes the values observed in flight (see table 3). The manufacturer s pred- 
ictions for the derivatives that are not functions of thrust coefficient are 
indicated by a single dashed line. These predictions are valid for both cruise 
and one-third flap down configurations. Details of the mathematics involved in 
the presentation of the manufacturer's derivatives are given in appendix B. 


LONGITUDINAL DERIVATIVES 


Figure 6 shows the MMLE-3 longitudinal derivative estimates for the Beech 
99 aircraft along with the corresponding predictions of the manufacturer. 
Comparison between the flight estimates and the manufacturers predictions 
shows close agreement except for the derivative . For some unexplained 

a 

reason, the flight data did not produce consistent estimates of this deriva- 
tive Other flight estimated derivatives are consistent, with the exception 
of the estimates of C and C resulting from one maneuver at an angle of 

q 6 e 

attack of 4.25 degrees. The plots of C N , C m , C N , and C are nearly 

a q 6 e 6 e 

constant with angle of attack, showing no observable effect of flap deflection 
All these parameters show good agreement with the manufacturer s estimates. 


LATERAL-DIRECTIONAL DERIVATIVES 


As with the longitudinal derivatives, the lateral-directional deriva 
tives are plotted with the manufacturer's estimates (fig. 7). The deriva- 


ti ves C V C V C„ p . C V C V C V C V y and are con- 

si stent with the derivative estimates of the manufacturer. C. is for 

6 a 


6 


the most part smaller than the manufacturer's estimates, while C 


n 


and C 


P 


V 


P 


are consistently larger in magnitude than the manufacturer s estimates for 
both cruise and one-third flap configurations. On the whole, the } a t®ral 
directional derivative estimates are repeatable and show consistent trends 
with angle of attack. 

CONCLUDING REMARKS 


A complete set of linear stability and control derivatives for a Beech 
99 airliner was determined using a modified maximum likelihood estimator. 

The derivatives were extracted for both the longitudinal and lateral -direc 
tional modes. The maneuvers were flown in smooth air at angles of attack 
ranging from 1.5 to 4.5 degrees. The first 29 maneuvers were flown in the 
cruise configuration and the last 27 were flown in a one-^hi rd ri 
configuration. The one-third flap down configuration had little effect on 
most of the stability and control derivative estimates. All 56 maneuvers 
produced satisfactory results. In general, derivative estimates from flight 
data for the Beech 99 airplane were quite consistent with the manufacturer s 
predictions, which are based on a combination of wind tunnel data, analytical 
estimates, and flight test data. 

Dryden Flight Research Center 

National Aeronautics and Space Administration 
Edwards, CA, November 27, 1978 
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APPENDIX A 


MAXIMUM LIKELIHOOD ESTIMATION PROGRAM AND 
EQUATIONS OF MOTION 

The analysis for this report was done with the MMLE-3 computer Drooram 
“ f the HMLE prograu (ref. 2). wlE- 3 is a gcncrTrnLiZme- 
lihood estimation program used at the Dryden Flight Research Center This 
section briefly describes the features of MMLE-3. 

feren«\ptwi!rti ike - lh0 °l eStiBWtes are deter ""'"ed by minimizing the dif- 
ference between the aircraft measured response and the calculated resDonse 

detemined by integrating the aircraft equations of notion. This rtlwTsl 

difference is formulated as a cost functional (discussed below) The .in- 

" ie Jrin t ’a,>r"« 0na ' I? * »0-ing the «ro~yna»U «?- 

ncienis in the aircraft equctions of motion. 

In general form, the equations of motion are: 

R(t) *(t) = A(t) x(t) + B(t) u(t) 
y(t) = C(t) x(t) + D(t) u(t) + Gq(t) 


tion| y of e o T r 0 Ce and R ; A H B, - ,D) ^ an b ® time functions because of the varia- 
.,I2 h io *5’ v ’ ?’ * d * dunng the maneuvers. Time varying matrices were not 

HgSre 5 a " alyS1S ° f th * BeeCh 99 data ’ exce P t the maneuve^shown in 
cost functional ,m likellhood estimates are obtained by minimizing the 


0 ( 4 ) 


= J jy c (t) - z(t)j ^GG*j *|yt) - z(t)J 


dt 


where 4 is the vector of unknowns, z is the measured response, and y f is the 

alQSrithmlrS?" 5 ^^ 5 ^ °J ** ™ LE T 3 uses a Newton-Balakrishnan iterative 
algorithm (ref. 2) to perform the minimization. 

cases ^averaoe 1 va US ? d in this rep0rt are given below - In man y 
matrices 1 ?K e ! I P aramete : s are u sed to obtain time invariant system 

but not for ofl Z USe sn,a11 ang1e a PP™ximations for p, 

Dut not for o0, or «j>. Symmetry about the XZ plane is assumed All enoies 

m these equations are in radians. The longitudinal state equations are 9 ' 


= - 9l 

mV 


[ 


(cos 0) (cos <>) (cos o) + (sin e) (si 


in a) 




/ 


a 


dl y » qSc C. 


0 = q(cos <j>) 


The longitudinal observation vector consists of the state vector concatenated 
with an observation of normal acceleration. The equation used for computing 
normal acceleration is: 


_ qS r A XAN . 

a = -S- C*. + q 

n mg N g M 


In expanded form, the longitudinal aerodynamic coefficients are: 


C N = C N a “ + C N 6 6 e + C Nq 
e 


C = C m a + c m §§ + C_ 6 + C 
m m a q 2V 6 e m 0 


C L = C N (for low a) 


The lateral-directional state equations are: 
p = jjj- C y + p(sin a) - r(sin a) + ^(cos 6) (sin 4») 

~ ^XZ = qsb C £ 

H X ' bI XZ = ^ sb C o 
$ = p + r(cos <J>) (tan 0) 

The lateral-directional observation vector consists of the state vector con- 
catenated with Tateral acceleration given by: 


* = as c ZAY w , 

y mg V g g 
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In expanded form, the lateral-directional aerodynami 


c Coefficients 


c v - y * Cy a ♦ cy r ♦ c Yo 

* c « r fs a a ♦ c « r » c 

a r u 

Cfl = S P * + C " r ^ * \ 6 a * V «r * C n, 


are 
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APPEND * a s* 


"MANUFACTURER'S DERIVATIVE ESTIMATES" 


Reference 5 contains the manufacturer's estimates of the total force and 
moment coefficients of the Beech 99 aircraft. These estimates are in an ana- 
lytical form which permits generation of stability and control derivatives by 
partial differentiation. These derivative estimates are valid for both the 
cruise and one-third flap configurations with the exception of the two esti- 
mates of C. . The resulting derivative estimates are functionally dependent 
a 

on the thrust coefficient. The derivative estimates and methods of the 
manufacturer are presented here to allow direct comparison of these estimates 
with flight results. 

The longitudinal derivative estimates are as follows: 

= 0.096 + 0.151 T c p^r deg (Cruise) 

- 0.101 + 0.151 T c per deg (1/3 flap) 

= 0.010 + 0.010 T c per deg 

= -0.034 - 0.032 T c per deg 

= -43.1 per rad 


3.40 per rad 


a 


a 


c/4 


m. 




The lateral-directional derivative estimates are as follows: 


0.014 - 0.0015 T c per deg 


-0.0023 per deg 

-0.010 per deg 

-0.0014 per deg 

0.00017 - 0.000022 a per deg 
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0.0025 per rad 



TOTAL THRUST 
where T = ; 

C qS 


-0.20 - 0.197a per rad 
+0.05 + 0.623a per rad 
0.394 per rad 
0.0079 - 0.762a per rad 
-0.515 per rad 
-0.199 - 0.385« per rad 
0.0027 per deg 
-0.0015 per deg 


At the steady state, low angle of attack conditions at which derivative 
extraction maneuvers were performed, the following approximations are valid. 



From reference 5: 
C D = 0.0275 

Thus , 

T = 0.0275 
c 


+ 0.0625 C L 2 


+ 0.0625 


(t-1 
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This equation shows that the cruise thrust coefficient is a function of ^ 

dynamic pressure, q, alone, because W and S are known constants. 


W = 39.59 kN 
S = 26 m 2 


'jbstituting these known values into the thrust coefficient equation, 

T, = 0.0275 ♦ 

C 0 2 


During the Beech 99 flight, dynamic pressure varied between 2 and 5 kN/m 2 . 
T c ranged from approximately 0 to 0.1. With this in mind, the Beech derivative 

estimates for the airplane were computed at the T c extrema, and these estimates 

were plotted with the flight- determined derivatives. 

Tht derivatives C and C must be resolved to a flight center of 

n R 

a c/4 P 

gravity position of 26 percent of the mean aerodynamic chord. This is ac- 
complished in the following equations: 


m 


= C 


m 

a 


(CG 


c/4 


flight 


- CG c/4 )' 


n„ 


= C 


P, 


Beech 


- C y (CG flight - CG c/4 )(c/b) 


For T c = 0, 

C ffl = -0.0289 per deg 

Of 

C n = 0.0014 per deg 


and for T = 0.1, 

= -0.0313 per deg 


C n = 0.00125 per deg 


^Reference center of gravity (CG^. . *) is 26 percent mean aerodynamic 
chord. ® 
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Wing: 


TABLE 1. GEOMETRIC CHARACTERISTICS OF 
BEECH 99 AIRPLANE 


Reference area, m 2 

Reference chord (mean aerodynamic chord), m 

Reference span, m 

Aspect ratio ..!!!,!*! 

Angl e of sweep , deg 

Ai ferons: 


26.01 

1.98 

13.98 

7.54 

0 


Area, (total )m 2 

Chord, m 

Span, (each) m 

Flaps: 


1.29 

0.28 

0.24 


Area, (total) m 2 
Span, (each) m . 
Chord, m 

Horizontal tail: 


Area, m 2 

Span, m 

Mean aerodynamic chord, m 

Aspect ratio 

Angle of sweep, deg 

Elevator: 


9.29 

6.82 

1.41 

5.0 

17.0 


Area, m 2 

Mean aerodynamic chord, m 

Vertical tail: 

Area, m 2 

Span, m ’ ' 

Mean aerodynamic chord, m 

Aspect ratio 

Angle of sweep, deg 

Rudder: 


4.17 

2.32 

1.92 

1.29 

19.5 


Area, m 2 

Mean aerodynamic chord, m 0 

Reference center of gravity, percent of reference chord ..26. 
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TABLE 3. INSTRUMENT LOCATIONS 


RELATIVE TO REFERENCE 


CENTER OF GRAVITY 



Instrument Distance Forward of Di stance below 

Reference CG, m reference CG, m 



a 


P 


7 

7 0 


a n 0. 381 

a x 0.102 

a y -1.5 0.102 


TABLE 4. INSTRUMENT RESOLUTIONS FOR BEECH 99 DATA 


Signal 

Resolution 

Of 

0.07 deg 

p 

0.08 deg 

p 

0.25 deg/sec 

q 

0.08 deg/sec 

r 

0.08 deg/sec 

6 

0.18 deg 

<t> 

0.35 deg 

a n 

O.Olg 

a x 

0 . 002g 

a y 

O.OOlg 

6 a 

0.08 deg 

6 e 

0.06 deg 

6 r 

0.10 deg 

V 

0.07 m/sec 


I 
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Case 

Dynamic 

Thrust 

Angle of 

Velocity, 

Flap 

Number 

Pressure, 

Coefficient 

Attack, 

m/sec 

Setting 


kN/m 2 


deg. 


deg. 

1 

4.73 

0.034 

1.7 

99 

0 

2 

4.62 

0.034 

1.7 

98 

0 

3 

4.56 

0.034 

1.7 

97 

0 

4/8 

4.77 

0.034 

1.6 

99 

0 

5/7 

4.68 

0.034 

1.6 

98 

0 

6/9 

4.67 

0.034 

1.6 

101 

0 

10 

4.95 

0.034 

1.9 

100 

0 

11 

4.56 

0.034 

1.8 

100 

0 

12 

4.61 

0.034 

1.8 

100 

0 

13 

4.63 

0.034 

1.8 

100 

0 

14/17 

4.50 

0.035 

1.8 

98 

0 

15/18 

4.69 

0.034 

1.8 

99 

0 

16/19 

4.53 

0.035 

1.8 

97 

0 

20 

4.83 

0.034 

1.5 

99 

0 

21 

4.75 

0.034 

1.6 

98 

0 

22 

4.67 

0.034 

1.5 

97 

0 

23 

4.65 

0.034 

1.5 

97 

0 

24/26/27 

4.80 

0.034 

1.5 

99 

0 

25/28/29 

4.80 

0.034 

1.5 

100 

0 

30 

2.23 

0.057 

3.2 

67 

15 

31 

2.20 

0.057 

3.4 

67 

15 

32 

2.20 

0.058 

3.3 

67 

15 

33 

2.17 

0.060 

3.5 

65 

15 

34/37 

2.10 

0.058 

3.3 

66 

15 

35/38 

2.25 

0.056 

3.1 

68 

15 

36/39 

2.22 

0.057 

3.2 

67 

15 

40 

2.04 

0.062 

4.2 

65 

15 

41 

2.07 

0.060 

4.4 

66 

15 

42 

2.09 

0.061 

4.4 

65 

15 

43/46 

2.08 

0.061 

4.3 

65 

15 

44/47 

2.03 

0.063 

4.4 

64 

15 

45/48 

2.04 

0.062 

4.3 

63 

15 

49 

2.24 

0.056 

3.0 

66 

15 

50 

2.20 

0.057 

3.1 

66 

15 

51 

2.15 

0.059 

3.1 

65 

15 

52/55 

2.14 

0.059 

3.2 

65 

15 

53/54/56 

2.16 

0.059 

3.2 

66 

15 
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TABLE 6. WEIGHTING VALUES USED IN ANALYSIS OF BEECH 99 DATA 


Signal 

Diagonal 

element of G INVERSE 

cr 

35 per deg 

q 

45 sec/deg 

e 

55 per deg 

a n 

300 per g 

P 

42 per deg 

P 

8.5 sec/deg 

r 

40 sec/deg 

♦ 

15 deg 

a y 

500 per g 















o 



(a) C N » 



Summary of longitudinal derivatives. 
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Figure 6 
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Figure 7. Summary of 1 ateral -directional derivatives 
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Figure 7. Continued. 
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